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Due to elastic anisotropy, two-dimensional patterning of substrates can promote weak azimuthal 
alignment of adjacent nematic liquid crystals. Here, we consider how such alignment can be achieved 
using a periodic square lattice of circular or elliptical motifs. In particular, we examine ways in which 
the lattice and motif can compete to favor differing orientations. Using Monte Carlo simulation 
and continuum elasticity we find, for circular motifs, an orientational transition depending on the 
coverage fraction. If the circles are generalised to ellipses, arbitrary control of the effective alignment 
direction and anchoring energy becomes achievable by appropriate tuning of the ellipse motif relative 
to the periodic lattice patterning. This has possible applications in both monostable and bi-stable 
liquid crystal device contexts. 


I. INTRODUCTION 

Surface anchoring, the promotion of a desired liquid 
crystal (LC) orientation at a surface [T], remains an im¬ 
portant problem in applications because precise tuning of 
anchoring parameters is often necessary for optimal de¬ 
vice performance [2j. Patterning the surface with a spa¬ 
tially varying preferred orientation is an attractive route 
to create alignment layers with desired anchoring prop¬ 
erties because both the effective anchoring strength and 
its orientation or easy axis can be altered by adjusting 
geometric features of the pattern P2 S]. Additionally, 
surfaces of appropriate symmetry [5] may promote mul¬ 
tiple stable easy axes leading to bistable devices HUE]. 
Bistability is desirable H0HI3] both for reduced power 
consumption and improved addressing of hi-resolution 
displays. Beyond displays, patterned LC systems are 
promising candidates as biosensors PUCE] and photonic 
devices pa na. 

Many methods exist to achieve patterning, encompass¬ 
ing both topographical and chemical approaches. These 
include mechanical rubbing m ng, photolithography 
mum], scribing with an atomic force microscope 
(afm) mam], microcontact printing of self-assembled 
monolayers (SAMs) [22H25] . and topographic surface fea¬ 
tures [26j. Since mechanical methods, such as rubbing, 
result in unwanted scratches or debris on the surface |27], 
and many methods do not scale well to high-volume man¬ 
ufacturing [25], SAMs have received much attention in re¬ 
cent years. Certain experimental methods show control 
over the azimuthal director angle as well as the polar 
anchoring angle mm - this is the focus of the work 
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presented in this paper. 

Striped surfaces, incorporating alternating regions pre¬ 
ferring planar and vertical alignment, have been well 
studied [f2M34] . For this pattern, the polar angle of the 
bulk LC is controlled by the average polar easy axis on 
the surface; the azimuthal alignment has an energy min¬ 
imum aligned parallel or perpendicular to the stripe ori¬ 
entation, depending on the ratio of the elastic constants 
[29]. Square checkerboard lattices are more complicated: 
an anchoring transition occurs in which the liquid crystal 
aligns with the lattice vectors for relatively strong surface 
anchoring, but switches to the diagonal for very weak 
anchoring [5j. Finally, both polar and azimuthal control 
over the bulk LC director orientation may be achieved 
with a rectangular checkerboard lattice. In this arrange¬ 
ment, certain rectangle ratios and anchoring strengths 
combine to shift the preferred director azimuthal angle 
from alignment with a rectangle edge to alignment diag¬ 
onally across the rectangle [35 j. 

For substrates constructed from squares or rectangles, 
the pattern element determines the symmetry and peri¬ 
odicity of the patterning. It is, however, straightforward 
to break this coupling by resorting to non-space-filling 
pattern basis motifs, such as circles or ellipses, and ar¬ 
ranging these on a periodic lattice. Importantly, this ap¬ 
proach provides a systematic approach by which to intro¬ 
duce the additional parameters needed to achieve truly 
independent control of polar and azimuthal anchoring an¬ 
gles and, so, set these at arbitrary target values. 

In this paper, we use a continuum approach to deter¬ 
mine the ground state of a liquid crystal in contact with 
a surface formed of a square lattice of circular or ellip¬ 
tical motifs. From this, we identify scenarios for which 
the basis symmetry and the lattice symmetry may fa¬ 
vor different alignments. This yields bistable configura¬ 
tions which are free from the constraints that pertain for 
square or rectangular motifs. The paper is organized as 
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Figure 1: Snapshots of Monte Carlo simulations of a nematic 
confined between periodically patterned circles, (a) Vertical 
circles on planar background for different values of thickness d 
and circle radius R. (b) Planar circles on vertical background 
for d/X = 1/2; the views show slices taken through the box 
mid-plane. 


follows: Monte Carlo simulation results with circle pat¬ 
terns are presented in section [TTJ In section pi} an analyt¬ 
ical continuum model is constructed for this arrangement 
with the simplifying assumption that the director lies at 
a constant azimuthal angle. We also construct a numeri¬ 
cal model that relaxes this assumption. Brief conclusions 


are presented in section IV 


II. SIMULATIONS 

The combination of Monte Carlo (MC) simulations 
and continuum theory has proven synergistic in previ¬ 
ous studies 0(35] because the alignment induced by a 
particular pattern depends dramatically on the length 
scales present: MC probes alignment around patterns at 
the order of a few molecular lengths while continuum 
theory permits modeling up to device scales. 

To gain an initial understanding of the effect of circle 
patterns on the adjacent nematic, we therefore first per¬ 
formed Monte Carlo simulations as described fully in [5]. 
Briefly, particle-particle interactions are modeled with 
the hard Gaussian overlap potential (HGO) and particle- 
wall interactions are modeled with the hard needle-wall 
(HNW) potential, where a hard axial needle of length 
ao k s is placed at the centers of particles of short diam¬ 


eter (To. By adjusting the needle length, both vertical 
(k s < 1) and planar alignment (k s > 2) can be achieved. 

Here, simulations of HGO particles of aspect ratio 
k = 3 were performed subject to confinement between 
two circle-patterned substrates separated by a distance 
d. Periodic boundaries were imposed in the x and y di¬ 
rections. Sharp boundaries were imposed between verti¬ 
cal and planar regions, needle lengths being specified in 
each region as described above. Particle configurations 
were initialized at low density and uniformly compressed 
in the x and y until orientational order was established. 
The corresponding in-plane box-length, A, then set the 
effective lattice periodicity of the chosen circular motif. 
An equilibration run of N = 500, 000 MC sweeps was 
conducted at each density, followed by a production run 
of a further N sweeps. 

Representative snapshots are displayed in Fig. [lja) 
for vertical circles on a planar background; correspond¬ 
ing plots for the reverse case are in Fig. 0b). For 
d = A = 4 Kao, the orientational ordering of the particles 
at the center of the film can be seen to depend on R, the 
radius of the circle. If d ~ A, an anchoring transition oc¬ 
curs with decreasing R. For large R, the film follows the 
alignment of the particles in the circle be that planar or 
vertical. However, a transition occurs with decreasing R, 
after which the orientation in the film becomes dictated 
by the pattern outside the circle. This behavior is ob¬ 
served for both vertical-on-planar and planar-on-vertical 
surfaces and essentially mimics that seen for other pat¬ 
terned films - the film orientation is dominated by the 
majority surface component. 

For thinner cells, where d = A/2; the nematic follows 
the pattern throughout the vertical distance for all values 
of R studied. This is similar to the “bridging” behavior 
observed in rectangular patterned systems [5j |35]. For 
the planar-on-vertical case [Fig. 0b)], the particle ori¬ 
entation on the top and bottom substrates at a vertical- 
planar boundary is particularly interesting, because the 
here the particles tend to align parallel to the boundary 
causing azimuthal distortion of the liquid crystal. 

A more subtle azimuthal transition is also apparent, 
however, from careful study of vertical on planar systems. 
As illustrated in Fig. [2j the preferred azimuth of the 
planar region is also dependent on R. Specifically, for 
larger circles particles in the in-plane region aligned along 
the box x and y- axes, whereas at smaller R they picked 
out the box diagonal. It is this azimuthal transition that 
we focus on in the remainder of this paper. 


III. CONTINUUM MODEL 

We now turn to the continuum analysis. The LC ori¬ 
entation is characterized by a director field, 


n = (cos 0 sin 0, cos 6 cos 0, sin 6 ), 


(i) 
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Figure 2: Apparent orientational anchoring transition, a) Top 
configuration for vertical on planar circles with R/X = 1/3 
and d/X = 1. b) Top configuration for R/X = 1/4. 


where 0 is the zenithal angle and is the azimuthal angle. 
The free energy of the LC will be given by the Frank free 
energy 


F = d 3 x [ffi(V • n) 2 + K 2 { n • V x n) 2 

+K 3 |n x V x n| 2 ] + J dSg{9 - 0 e ), (2) 

s 



Figure 3: (a) Schematic of the unit cell domain with im¬ 
portant length scales labeled for a circle-patterned surface 
(a = b = R and uj — 0). (b) Definition of the polar angle 6 
and the azimuthal angle </> of the LC director. 


and we introduced the dimensionless polar anchoring pa¬ 
rameter 


with a harmonic anchoring potential 


g(o-o e ) = ^(o\ z=±d/2 ~e e ) 2 . (3) 

The coordinates are set up as follows: consider a unit 
cell defined on the box with corners at (0,0, —d/2) to 
(A, A,+d/2). Each surface at ±d /2 contains an ellipse 
centered on the unit cell with semi-major axis a oriented 
at an angle uj with respect to the x-axis and semi-minor 
axis b. The surfaces promote homeotropic (6 = 7t/2) 
alignment within the ellipse and planar alignment (0 = 0) 
outside, as shown in Fig. [3] 


A. Solution 


Le 


K i 
WqX' 


(5) 


Within this approximation, minimization of the free 
energy yields an anisotropic Laplace equation, 


V T • A • V0(x, y , z) = 0, (6) 

where a = 1 + r and /3 = 1 — r . This may be solved 
using the series, 


oo 

9(x,y,z)= 2A nm cosh(r nm z)e i2w ( nx+my \ (7) 

n,m= — oo 


Following previous work, we make the two-constant 
approximation by setting K\ = K% and K^/K 1 = r |29]. 
Additionally, if the polar angle depends on all of the coor¬ 
dinates 6 = 0(x,y, z) and the azimuthal angle is constant 
in space (/> = (/>o [35], the bulk free energy may be rewrit¬ 
ten as a quadratic form, 

f =*yJ ^0f-MV0)d 3x +d/±^ I dS(0\ z=±d/2 -e e ) 2 

S 

(4) 

where 

— cos20 o ) ^/3sin20 o 0\ 

A — [ |/3sin20 o | (ck + /3 cos 20o) 0 ] 

V 0 0 1/ 


where 

Vn / m = \Ja(m 2 + n 2 ) + f3(m 2 — n 2 ) cos 20 O + 2/3nm sin 20 o . 
\/27r 

( 8 ) 

In order to satisfy the boundary conditions, the patterned 
easy axis is first expanded in a Fourier series, 


0e(x,y)= Snme i2 ^ nx+my) . (9) 

n,m=0 

To determine the coefficients S nm , we assume that the 
background surface promotes 6 e = 0 while the elliptical 
pattern promotes 0 e = 7r/2, henceforth referred to as 
“vertical on planar” patterning; the alternative “planar 



















on vertical” arrangement is trivially obtained from this 
solution by making the substitution, 
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6 7r/2 — 0, (10) 

which leaves the energy invariant. The S nm are therefore 
evaluated by integrating, 

$nm = 7^ I exp (27r i(nx + my))dxdy , (11) 


0 C = 


oz 


z=±d/2 


(14) 


Here, the ± refers to the direction of the outward normal 
to the LC boundary. Inserting Eqs. 0 and 0 into ( pH| , 
we obtain, 


over a domain V defined by the ellipse equation, 


(x - x c ) T R(u) ^ 1// Q a y h 2 ^ R(-u) (x - x c ) < 1. 

( 12 ) 

Here, R(u) is the 2D rotation matrix and x c = (|, \) is 
the center of the ellipse. As shown in the appendix, the 
integral © can be performed analytically to yield, 


c _ dbir , x n _|_ m 

dnm — 2 V D 


•A (2?r y a(j 2 m + 6; 2 m ) 




2 

nm 


b' 2 

nm 


(13) 


A rt.m, — 


2 (Ler nm sinh (^y) + sinh(%y)) 


(15) 


Note that as —>> 0, this recovers the rigid anchoring 

condition. The solution for Q(x,y,z) is now obtained by 
inserting Eqs. (15), (13), and (J8| into Eq. 


As for other patterns 0|35], the director follows the 
surface pattern at ±d/ 2, while relaxing to a uniform ori¬ 
entation, equal to the average polar angle promoted by 
the surface far away from the boundaries. 


where J\{x) is a Bessel function of the first kind, a nrn = 
(n cos uj + m sin uj)a, and b nm = (■ n sin uj — m cos uS)b. 

Having expanded the easy axis in a suitable form, 
the coefficients A nm can be determined by imposing the 
Robin boundary condition at did/2. 


B. Circular patterns 

For a circular surface pattern, set a — b in Eq. 
Evaluation of the volume integral in Eq. 0 then yields 
an expression for the bulk energy of the LC, 


F b = 


oo 

E 

n,m= —oo 


2KiA 


2 

nm 


B 1 A 


2 diag(l, 1, —1) + 


sinh ( Tu rnd) 
2r nn 


h)B 


(16) 


where 

/ 27m \ 

B = I 2irm J . 

y ^nm J 

Similarly, evaluation of the surface integral over the sur¬ 
faces at +d/2 and —d/2 gives the surface energy of the 
LC, 


_ A'i(S„ m -2yl n „co8h(^)) 2 

F ’= ^ - T, -' 

n,m= —oo 

The bulk and surface energy of the LC are shown in 
Fig. [4] as a function of the azimuthal angle <p for fixed cir¬ 
cle radius and polar anchoring energy. The bulk energy 
always prefers alignment along the x- or y- axis, but with 
decreasing strength as the unit cell thickness decreases. 


Meanwhile, the surface energy prefers director alignment 
at a 45-degree angle to the axes. This surface preference 
grows slightly stronger as cell thickness decreases. Re¬ 
sults for planar on vertical patterning are identical due 
to the invariance of the energy under the linear transfor¬ 
mation (10). 


We estimate W e ff , the effective azimuthal surface an¬ 
choring energy m as the energy difference between the 
0o = 7 t/ 4 and the 0 0 = 0 states per unit area. A positive 
value of W e // indicates a preference for alignment along 
the x- or y- axis, while a negative value of W e ff indicates 
a preference for diagonal alignment. Fig. |5[a) shows 
the effective surface anchoring potential as a function of 
We, the polar anchoring energy, for a series of cell thick¬ 
nesses. The surface contribution dominates for a thin cell 
with strong anchoring Wq ~ 100, which corresponds to 
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Azimuthal angle 0 O 


Figure 4: (a) LC surface energy (dotted), bulk energy 

(dashed), and total energy (solid) as a function of the az¬ 
imuthal angle for a unit cell of thickness d/X = 1.0 (brown) 
and d/X = 0.25 (cyan) for planar on vertical patterning, 
r /X — 0.5, and L = 0.01. For the thicker unit cell, the surface 
energy azimuthal preference is much weaker than that of the 
bulk and the overall preference aligns with that of the bulk. 
However, thinner cells exhibit weaker bulk energy preference 
and the total energy favors alignment along 0o = 7 t/ 4. (b) 
Corresponding calculated structures for i) d/X = 1.0 and ii) 
d/X = 0.25. 


Lq = 0.01. In the limit of rigid anchoring (Wq oo) 
or extremely weak anchoring (Wq —»• 0), the effective an¬ 
choring vanishes. For weak anchoring, the nematic effec¬ 
tively ignores the pattern, while for rigid anchoring the 
surface follows the prescribed pattern exactly. 

To determine the parameter space in which parallel 
and diagonal alignment are each favored, we display in 
Fig. |5jb) the strength of the azimuthal energy preference 
computed from the combined bulk and surface energies, 
shown for Lq = 0.100, and Lq = 0.001. The diagrams 
indicate two regions in which the the surface preference 
overrides the bulk preference for a sufficiently thin cell. 
The strongest of these regions is the thin brown section 
at R/X = 0.5, the largest possible circle radius. This 
region of parameter space is dominated by the surface 
term, likely due to the inability of vertically/horizontally 
aligned LCs to effectively fill space between two abutting 
circles whose edges are nearly perpendicular to the LC 
director. For circles much smaller than R = 0.5, the 
preference becomes very weak indeed. 

The comparable magnitude and conflicting preferences 
of the bulk and surface energies at these length scales 
suggests that the constant-0 approximation may be too 
restrictive for this system. Unlike the rectangular and 
square patterns previously considered, the alignment di¬ 
rection is promoted exclusively by the lattice while the 
basis favors no particular alignment. The results of the 
Monte Carlo simulations also suggest this: in the configu¬ 
rations shown in Fig. 0b) , the particles tend to align tan¬ 
gentially around the edge of the circle because for r < 1 
the energetically cheapest way to achieve the vertical-to- 




Figure 5: (a) Effective anchoring energy as a function of polar 
anchoring energy Wq for R/X — 0.5 and several cell depths. 
As the cell depth decreases, the surface preference becomes 
more pronounced for a given value of Wq. (b) Phase dia¬ 
gram showing the strength and orientation of the preferred 
azimuthal alignment angle as a function of cell depth d/X 
and circle radius R/X for i) Lq — 0.100, and ii) Lq = 0.001. 
Brown regions indicate diagonal alignment 0o = 7r/4 while 
cyan regions regions prefer alignment with the lattice vectors, 
0o — 0 or 0o = 7 t/ 2. The brown regions expands into larger 
cell depths as the anchoring strength increases. 


planar transition around the perimeter of the circle is 
through a twist deformation. In section liilDl therefore, 
we numerically minimize the free energy (|2| with respect 
to a completely arbitrary director profile to quantify the 
effect of azimuthal variations. 


C. Elliptical patterns 

We now consider elliptical patterns. For long ellipses, 
one might expect the effective azimuthal alignment to lie 
parallel to the semi-major axis, resembling the situation 
with alignment on striped surfaces|29j. Hence, by ad¬ 
justing 0, it should be possible to control the preferred 
azimuthal angle arbitrarily and, by tuning the aspect ra¬ 
tio, also control the effective azimuthal anchoring energy. 
The control parameter space to consider is greatly ex¬ 
panded: while the cell depth d/X and the anchoring po¬ 
tential Wq remain parameters, the circle radius R/X is 
replaced by the semi major axis length 6/A, aspect ratio 
6/a, and the alignment angle uj. From the structure of 
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Azimuthal angle 0o 



Figure 6: Total LC energy as a function of the azimuthal angle 
shown for (a) several ellipse rotation angles with b/a = 1.2 
and (b) for several different aspect ratios with uj = 7t/4. Both 
panels have d/X = 1, Lq — 0.01, and a coverage fraction 
equal to that of a circle with R/X = 0.4. Rotation of the 
ellipse results in migration of the preferred azimuthal angle 
from alignment with the axes of the unit cell to alignment 
with the semi major axis of the ellipse. However, panel (b) 
shows that this transition in azimuthal angle preference is not 
immediate, but instead passes through several smaller angles 
as aspect ratio increases. 

the solution 0, we expect the cell depth and anchor¬ 
ing potential to have similar effects in both patterns and 
so we focus on the effects of the new parameters in this 
section. 

In Fig. [6j we show effective azimuthal anchoring pro¬ 
files W e ff((f)) calculated for a variety of values of uj and 
b/a. The coverage fraction of the pattern, or equivalently 
the area of the elliptical motif, is kept constant in order 
to fix the effective polar angle. An immediately obvious 
feature, compared to the equivalent profiles for circular 
patterns plotted in Fig. [4j is that the mirror symmetry 
of the pattern about </> = 7r/2 is entirely broken, leaving 
behind a non-symmetric anchoring potential reminiscent 
of the structures fabricated in m- 


For fairly modest aspect ratios, the alignment angle 
of the ellipse controls the preferred azimuthal angle such 
that the energy minima occur at 0o = w. For smaller 
values of a;, or aspect ratios close to unity, the results 
are more complex: For instance, the energy minimum for 
uj = 7r/8 in Fig. [6^ a) occurs at 0o = 7r/16 instead of <j )o = 
7r/8. Also, in Fig. mb) we see that, for a rotation angle 
of 7t/4, the azimuthal preference moves from alignment 
with the sides of the unit cell to alignment with the semi 
major axis of the ellipse gradually, preferring </>o = n/16 
for an aspect ratio of 1.05 and 0o = tt/8 for an aspect 
ratio of 1.1. 

The mechanism for this transition is two-fold. Firstly, 
the surface energy consistently prefers azimuthal align¬ 
ment along the semi-major axis of the ellipse. Though 
the angle preferred by the surface remains the same for 
any non-unit aspect ratio, the strength of the preference 
grows with increasing aspect ratio. Second, as the as¬ 
pect ratio grows, the bulk energy of the LC tends to pre¬ 
fer an azimuthal angle aligned with the ellipse alignment 
angle, but this move happens slowly such that small-to- 
moderate aspect ratios result in a preferred angle some¬ 
where between 0o = 0 and 0o = c u. The magnitude of 
the bulk energy preference also grows with increasing as¬ 
pect ratio, but less dramatically than that of the surface 
energy preference. 

The preferred angle 0o is therefore the result of a ten¬ 
sion between surface and bulk effects. We display the 
direction and strength of the effective anchoring poten¬ 
tial as a function of aspect ratio in Fig. 0a). Fig. 0b) 
presents a phase diagram for preferred angle as a function 
of the semi-major axis and aspect ratio of the elliptical 
motifs while holding ellipse rotation angle, cell depth, 
and anchoring potential constant. As expected, increas¬ 
ing aspect ratio leads to alignment of the effective easy 
axis with the ellipse long axis. Less intuitive is the fact 
that smaller ellipses transition from alignment with the 
lattice to alignment with the basis at lower aspect ratios. 
The decreased scale of the bulk LC energy for small sur¬ 
face features drives this as unit cells with small motifs 
have surface and bulk energies at the same scale. Thus, 
increases in aspect ratio quickly shift the LC azimuthal 
preference to align with the basis, overcoming the influ¬ 
ence of the bulk. 


D. Numerical model 

As briefly discussed earlier, the Monte Carlo simu¬ 
lations produced configurations in which the particles 
tended to align tangentially with the boundary of the 
circle pattern. These are shown in Fig. [8] and, together 
with the results of the previous sections, suggest that the 
constant 0 approximation may be too restrictive for this 
system. We therefore performed a numerical minimiza¬ 
tion of the Frank energy 0 again using the two-constant 
approximation but now allowing the director to vary ar¬ 
bitrarily in 3D. A Cartesian representation of the director 
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a) b) 




Figure 7: (a) Magnitude (dashed, left axis) and direction (solid, right axis) of azimuthal angle preference given a growing aspect 
ratio for ellipse alignment angles of u — 7t/4 (cyan) and uj — 7r/8 (brown) with d/X = 1, Lq — 0.01, and a coverage fraction 
equal to that of a circle with R/X = 0.4. (b) Phase diagram of azimuthal angle preference of the entire unit cell as a function 
of aspect ratio b/a and semi major axis length b/X with uj = 7r/4, d/X — 1, and Lq — 0.01. Note these parameters do not 
maintain a constant coverage fraction. 



Figure 8: (a) Results from 3D director minimization model. 
Shown is the director orientation on the bottom face of the 
unit cell for planar circle patterns of radius R/X = 0.5, R/X = 
0.4, R/X — 0.33, and R/X — 0.2. These simulations are 
performed with Lq = .001, d — 1.0, a final grid size of 20 x 
20 x 21, and an initial guess that had the director aligned 
along the y- axis, (b) Corresponding results from Monte Carlo 
simulations. 


n = (n X: n y ,n z ) was used and the unit length constraint 
h • n = 1 enforced locally. The energy was discretized 
using second-order finite differences and minimized using 
an adaptive gradient-descent relaxation method with line 
searches. To improve convergence, successive refinement 
was used: an initial guess is relaxed on a coarse grid, then 
interpolated and relaxed onto successively finer grids. At 
each step, the system was relaxed until the energy con¬ 
verged. 

Results from the relaxation model for a set of planar 
patterned surfaces are shown in Fig. J^a). An initial 
guess with the director aligned along the y -axis was used. 
After relaxation, the director adopts an orientation tan¬ 
gent to the pattern edges, as seen in the correspond¬ 


ing Monte Carlo simulations, Fig. [8^b). The tangen¬ 
tial alignment arises because it corresponds to a twist 
deformation across the vertical planar boundary, which 
is energetically cheaper than a bend or splay deforma¬ 
tion. The behavior breaks down for smaller circle radii 
because the bend deformation required to follow the arc 
of the feature becomes too energetically expensive. 


IV. CONCLUSION 

This paper has considered the alignment behavior of 
square arrays decorated with elliptical motifs, demon¬ 
strating that such patterns can be used to create surfaces 
with controllable anchoring energy and easy axis by vary¬ 
ing the period of the pattern, ellipse orientation, aspect 
ratio, and coverage fraction. Given the two-constant ap¬ 
proximation and the assumption of a constant azimuthal 
angle, the director configuration and energy can be com¬ 
puted analytically. Depending on the circle radius, cell 
depth, and polar anchoring strength of the alignment ma¬ 
terial, the ground state may have azimuthal alignment 
along either of the lattice vectors, or diagonally. 

Our study offers invaluable advice for applications be¬ 
cause, while the elliptical patterned surface offers a re¬ 
markable degree of control over the anchoring proper¬ 
ties, the design parameter space for the pattern is large. 
Briefly, surfaces patterned by rotated ellipses allow con¬ 
trol of the azimuthal angle over a continuum of values 
between the lattice vectors and the diagonal depending 
on the orientation of the ellipse. In these cases, the ellipse 
aspect ratio controls the effective anchoring energy. The 
behavior is nontrivial, however, and includes regions of 
bistability, as well as an azimuthal anchoring transition 
for some designs. 




























































Unlike for other patterned surfaces previously studied, 
the constant azimuthal angle approximation is of only 
limited use. A numerical study showed a preference for 
tangential alignment along the vertical-planar boundary 
of the pattern in excellent agreement with Monte Carlo 
simulations. 

While the present study has considered a flat surface, 
the results may also be important for experimentalists 
using arrays of posts to align liquid crystals as has been 
reported in da OH ED- It seems that only a very mod¬ 
est amount of anisotropy in the shape of the pillars, for 
example using posts of elliptical cross section, may break 
the square symmetry and lead to significantly better az¬ 


imuthal alignment if desired. 
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Appendix 


To evaluate the Fourier coefficients for the easy axis, first substitute x 
and p| ) to be rewritten as 


x' + 1/2 and y y' + 1/2, which allows 


Sum = ^ exp (7 ri(n + m)) JJ exp ( 27 t i(nx' + my'))dx'dy ', 


(17) 


and 


/ I / • \2 / / • / \2 

x cos uj + y sin uj \ ( x sin uj — y cos uj 


< 1, 


(18) 


respectively. Note that the exponential pre-factor in ( |17| ) is (—l) n+m since n and m are integers. 

Next, rotate the coordinates via the transformations (x'coscj + y' sincj)/a x" and (x'sincj — y r cosuS)/b y ", 

and integrate using these new coordinates, 


Sr>/m — 


abir 


(—l) n+m JJ exp ^27ri((ncoscc;+ msinuj)ax" + (n since; — m cos cj) by ")^ dx"dy". (19) 


(x")2 + (i/ // ) 2< 1 


Define a' = (ncosu; + msinu;)a, b' = (n since; — mcosu;)6, a = (a', 6'), and g(s) = exp(is) and convert to polar 
coordinates such that (x",y") —>> rf where £ = (cos#, sin#): 


Sn/rn — 


abir 


(—1 ) n + m y" r dr J g(2itra • f)df. 


Eqs. (1.2) and (1.5) from [40] allow us to evaluate the inner integral to find 


Snm — 


abir 


r ^ 

(—l) ra+m J 2irrJo ^27r r\J a' 2 + 6' 2 ^) dr. 


Finally, let k = 2ir\/ a' 2 + 6' 2 and s = kr, and note that f xJo(x)dx = xJ\(x) to obtain, 


( 20 ) 


( 21 ) 


Sum = a ^(-l) n+m ^J\j 0 (s)ds 

= ^(-l) n+m yJi(*0 

afoy n n+m J i (Wfl /2 + & /2 ) 
2 1 J Va' 2 + &' 2 


( 22 ) 
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the result stated in Eq. (13). 
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